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ABSTRACT 

In recent years, upper limits on the cosmic microwave background (CMB) anisotropies combined with 
predictions made by theories of galaxy formation, have been extremely powerful in ruling out purely 
baryonic dark matter (BDM) universes. However, it has recently been argued that the absence of a 
prominent second peak in the anisotropy spectrum measured by the BOOMERanG-98 and MAXIMA- 1 
experiments seems to favour a ABDM model when combined with standard Big Bang Nucleosynthesis 
(BBN) constraints. In this Letter, we further investigate this result showing that, using the CMB data 
alone, a purely baryonic adiabatic model of structure formation seems unlikely if the universe is flat 
(CI = 1). Combining the CMB data with supernova type la (SNIa) data renders purely baryonic models 
inconsistent with flatness at high significance and more than 3cr away from both the BBN constraints 
and the HST key project result h = 0.72 ± 0.08. These results indicate that only a radical revision 
of cosmology with ad hoc properties could bring baryonic models such as those advocated by MOND 
enthusiasts back into agreement with current observations. 

Subject headings: cosmology: cosmic microwave background — cosmology: dark matter 



1. INTRODUCTION 

Since the pioneering work of Zwicky (1933), astronomers 
have found compelling evidence to suggest that the ma- 
jor contribution to the overall mass density of the uni- 
verse is in the form of dark matter. Even if its p recise 



reasoning is adopted, currently discussed versions of the 
present cosmological scenario, incorporating a fine-tuned 
cosmological constant and dark matter particles with ad 
hoc properties, lose the comp elling aesthetic simplicity of 



nature remains unknown (see e.g. Bergstrom 200C for a 



the original CDM model (see |5ellwood fc Kosowsky 2000 
for a recent review). Of course , other, more astrophysi 



recent review), the currently favoured hypothcoio of cold 



cal, r esolutions may be possible Binney, Gerhard and Silk 



dark matter (CDM) is in good agreement with an onor - 

mous body of data: anisotropies in the cosmic microwave 
background (CMB), large-scale structure surveys, cluster 
abundances, structure of the L yman-q forest and a host 
of other measurements (see e.g. Bahcall et al. 1999). Fur- 
thermore, in the past 30 years, upper limits on the CMB 
anisotropies and predictions made by theories of galaxy 
formation have been extremely effective at ruling out mod- 
els based on purely baryonic or hot dark matter (see e.g 



2001 



Doroshkcvich, Zel'dovich fc Sunyaev 1978 ; Wilson fc Silk 



1981|; |Kaiser 1983fc |Whitc, Frcnk fc Davis 1983j J 

The quality of astrophysicai data has improved and it 
currently seems likely that the present fractional overall 
energy density fl includes contributions from a cosmologi- 
cal constant A (Qa) and from cold (fi c dm), hot (Clhdm) and 
baryonic (Clt) dark matter. Also, it appears that the den- 
iitioD of thooo compononto are all within one or two ordcro 



of magnitude of each other, — While the similar densities 



The recent CMB anisotropy measurements from the 
BOOMERan G-98 (|dc Bcrnardis et al. 2000| ) and 
MAXIMA-1 (Hanany et al. 2000) experiments have pro- 
vided new insights into the cosmological parameters. For 
the first time, the anisotropy angular power spectrum has 
been measured over a wide range of angular scales, from 
multipole I ~ 50 up to I ~ 800 with errors of the order 
of 10%. The data sets confirm that there is a peak in the 
angular power spectrum at I ~ 200 with a steep decline 
in power from I ~ 200 to I ~ 300. While the presence 
of such a peak is consistent with the predicted acoustic 
oscillations in the adiabatic inflationary scenario, the ab- 
sence of prominent secondary peaks after £ ~ 300 is an 
unexpected result, suggesting that the value for the phys- 
ical baryonic density is ~ 50% higher than that predicted 
by Big Bang Nucle o synthesis (BBN) (see e.g. Tegmark 
& Zaldarriaga 2000J paffe et al. 2000| ). We note however 



of dark 'energy' and 'matter' can possibly be explained 
by advocating a q uintessence scalar field which 'tr acks' 
the matter density (Steinhardt, Wang & Zlatev 1999), the 
reason why the densities of the 3 matter components are 
so similar is not so apparent. 

Predictions made by structure evolution theories based 
on the ACDM model disagree with some galaxy obser- 
vations, which suggests that, for the model to hold, as- 
sumptions made about the properties of CDM particles 
must be modified ( see e.g. [Bpcrgcl fc Steinhardt 1999 



that this discrepancy with the standard model is no more 
than 2 a w hen the most recen t CMB anisotropy results 
arc included Mason et al. 2001. 



Burkert 200C ; Madsen 2000; |Moore et al. 2000). If this 



In recent months, several solutions have been advocated 
to reconcile BBN with the larger values of the baryon den- 
sity suggested by the CMB data. In a recent letter by 
McGaugh (2000), the simple and controversial solution of 
considering only baryonic dark matter (BDM) plus a cos- 
mological constant has been proposed. The lower CMB 
second peak can thereby be explained by the collisional 
damping expected in a BDM universe, while maintaining 
the baryon density at a level compatible with the BBN 
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constraint. As is well known, purely baryonic adiabatic 
models are in severe contradiction with the measured am- 
plitude and shape of galaxy and cluster correlations. How- 
ever, there remain serious issues about the possibly scale- 
dependent bias parameter. Hence McGaugh's result is in- 
teresting for various reasons: it resolves the disagreement 
with BBN by using a simpler model (fl c dm — ^hdm = 0), 
it demonstrates the model-dependence of theoretical fits to 
the CMB data by exploring a previously uncovered corner 
of parameter space and, finally, it questions our theoretical 
understanding of the structure formation process. 

If it turns out that the CDM hypothesis does not work 
well at some cosmological epochs and scales, then even the 
exotic scenario of a purely baryonic model with modifica- 
tions to standard gravity could become competitive, espe- 
cially in view of the insistence of MOND proponents that 
this alternative theory of gravity merits serious examina- 
tion ( [McGaugh fc do Block 199l| ganders 2000| ). Thus, it 
is timely to study how compatible the CMB data is with a 
BDM adiabatic universe, what consequences this assump- 
tion can have on the remaining cosmological parameters, 
and the need for CDM from independent observations of 
galaxy clustering. 

In this Letter, we perform an analysis in cosmological 
parameter space. Using the recent CMB data alone we 
show that a purely baryonic adiabatic model appears un- 
likely if the universe is flat (f2 = 1). By further combining 
these results with supernova type la (SNIa) data, we rule 
out baryonic adiabatic models with high confidence. 

2. METHOD 

The structure of the Cg spectrum depends es sen- 



Sugiyama & Silk 199?; 


Efstathiou & Bond 1999 


chiorri k Griffiths 2000 


and references therein) 
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physical baryonic density u>i, = flbh and the overall 



matter density uj m — Q m h = (f2. 



cdm + ttb)h 2 which 
define the size of the acoustic horizon at decoupling, 
and the 'shift' parameter R related to the geometry 
of the universe through R — 2y/\£lk\/£l m /x(.y)i with 

y = vW/ Zd "p m (i + ^) 3 + ^(i + ^) 2 + ^A]" 1/2 ^, 

where the function x(v) depends on the curvature of the 
universe and is y, sin(y) or sinh(y) for flat, closed and open 
models respectively. 

A decrease in w cc j m , with u>b and R remaining con- 
stant, will reduce u> m . This decreases the redshift 
of equality causing the peak in the spectrum to be 
shifted towards smaller angular scales. Constant u> m 
can be maintained by increasing cj&. However, an in- 
crease in lui, will decrease the sound speed at decoupling 

--dec ))) again reducing the 
acoustic horizon size and shifting the peak to higher £'s. 

To keep the position of the first peak fixed while u>b and 
uj m are varied, R has to be appropriately tuned. This 
can be achieved by increasing f^A- The shift parameter 
R is sensitive to variations in Q\ as Vt m — > 0, and for 
fi TO = S^fc ~ 0.05 a 10% increase in fl\ can produce a 50% 
decrease in £ pea k- Thus, a purely baryonic model can pro- 
duce the observed peak structure at £ ~ 200 only if R is in- 
creased with respect to the corresponding flat CDM model. 
We can therefore anticipate that viable purely baryonic 
models will be A-dominated with closed geometries. 



Motivated by these considerations, we compare recent 
CMB and SNIa observations with a set of models with 
parameters sampled as follows: Q m = Qt, — 0.015, ...,0.5; 
0\ = 0.80, ...,1.04 and h = 0.40, 0.95. We vary the 
spectral index of the primordial density perturbations 
within the range n s — 0.50, 1.50 and we rescale the fluc- 
tuation amplitude by a pre-factor Cio , in units of Ci(f BE . 
The theoretical models are c omputed using the publicly 
available CMBFAST program (Seljak & Zaldarriaga 1996) 
and are compared with the recent BOOMERanG-98 and 
MAXIMA-1 results. The power spectra from these ex- 
periments were estimated in 12 and 10 bins respectively, 
spanning the range 25 < £ < 785. In each bin, the spec- 
trum is assigned a flat shape, £{£ + l)Cg/2ir = Cb- 

Following Bond, Jaffe & Knox (1999) we approximate 
the signal Cb inside the bin to be an offset lognormal 
distribution, such that the quantity Db — ln(Cs + xb) 
(where xb is the offset correction) is a Gaussian variable. 
The likelihood for a given cosmological model is then de- 
fined by -21ni = (L>f - D^)M B b' (Dg, - £>ff,) where 
Cg (C B X ) is the theoretical (experimental) band power, 
Xb is the offset correction and Mbb' is the Gaussian cur- 
vature of the likelihood matrix at the peak. We consider 
10% and 4% Gaussian distributed calibration errors for 
the BOOMERanG-98 and MAXIMA-1 experiments re- 
spectively. We also include the COBE data using Lloyd 
Knox's RADPack packages. Proceeding as in Dodelson & 
Knox (2000), we attribute a likelihood to a point in the 
(Q m = fib, ^a) and (fib, h) planes by finding the remain- 
ing 3 parameters that maximise it. We then define our 
68%, 95% and 99% contours to be where the likelihood 
falls to 0.32, 0.05 and 0.01 of its peak value, as would be 
the case for a two dimensional multivariate Gaussian. 



3. RESULTS 

In the top panel of Figure 1, we plot likelihood contours 
in the (f2 m = £Ia) plane by applying the maximiza- 
tion/marginalization algorithm described above and using 
only the BOOMERanG-98 data and limiting our analy- 
sis to models with age to > lOGyr. As expected, the low 
sound speed at decoupling due to the high baryon content 
makes only closed models compatible with the observa- 
tions. The deviation from flatness becomes less and less 
important as Qb — > 0, as one would expect from the ex- 
pression for R. However, the decrease in the redshift at 
equality causes positive curvature models to be preferred. 

In the parameter range we are sampling, we find that 
flat models are excluded at ~ 2a. In order to further test 
the discrepancy with a flat universe we also include CDM 
models in the analysis with f2 c dm = 0.05, 0.45. Compar- 
ing our models to the BOOMERanG-98 data alone and re- 
stricting the analysis to fl = 1, we find that h 2 il c d m > 0.04 
at the 95% C.L. Thus flat BDM models are excluded with 
significance. 

For the BOOMERanG-98 analysis with tt cdm = 0, the 
best-fitting model is a closed model with Qb = 0.135, 
n A = 1.0, n s = 0.94, h = 0.45, C w = 0.5 and the 
BOOMERanG-98 calibration left untouched. This pro- 
vides a good fit to the data (see Figure 2, top panel). In- 
cluding the MAXIMA-1 data points produces a different 
best-fitting model: Q, B = 0.035, fl A = 1.00, n s = 0.92, 
h = 0.85, Cio = 0.5 with a ~ 12% upward calibration 
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FlG. 1. — The likelihood contours in the (£Ib, Ha) plane, with 
the remaining parameters taking their best-fitting values. Only 
the BOOMERanG-98 data is considered here. The contours 
correspond to 0.32, 0.05 and 0.01 of the peak value of the like- 
lihood, which are the 68%, 95% and 99% confidence levels re- 
spectively. The dashed line going from left to right diagonally 
across the plot indicates flat models. The top panel is with a 
prior on the age of the universe t > 10 Gyr. The bottom panel 
includes the BBN prior n b h 2 = 0.019 ± 0.002 and the SNIa 
data. The results of these two priors are incompatible at ~ 3a. 



for the BOOMERanG-98 points and a ~ 5% downward 
calibration for MAXIMA- 1. 

As the low value of C%o already suggests, these best- 
fitting models are in contradiction with the COBE data. 
Including COBE in the BOOMERanG-98 analysis does 
not significantly affect the best-fitting cosmological param- 
eters which remain tts = 0.135, fl\ = 1.0, ns = 0.94, 
h = 0.45, but the model now requires Cio = 0.8 and 
an upward shift in the calibration for BOOMERanG-98 
of - 30%. This model is ruled out by the COBE and 
BOOMERanG-98 data at the 95% C.L., although this dis- 
crepancy can be resolved by adding a gravitational wave 
(GW) contribution on large scales. The top panel of Fig- 
ure 2 indicates that the inclusion of a GW component such 
that CJq — Cf leaves the small scale behaviour practi- 
cally unchanged but provides enough CMB power on large 
scales to match the COBE normalization. 

A GW contribution as large as Cj/Q? ~ 1 is incom- 



patible at the 95% C.L. with flat ACDM models when 



the B BN prior is assu med in CMB analyses (Kinney, Mel- 
chiorri & Riotto 2000). GWs leave a characteristic imprint 
on polarization power spectra, via fluctuations in the B 
magn etic-type-parities channel ( Kamionkowski, Kosowsky 
& Stebbins 1996, Seljak & Zaldarriaga 199C) that vanish 
in the case of scalar fluctuations. However, the expected 
polarization amplitude for BDM models is very low, below 
0.5/iif, and future experiments will not have the sensitivity 
to detect the B channel signal unless a substantial amount 
of reionization occurred in the past flScljak fc Zaldarriaga 
19960 
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FIG. 2.— In the top panel, the BOOMERanG-98 and the 
COBE data are plotted together with the best-fitting baryonic 
model to the BOOMERanG-98 data. In order to fit the COBE 
data, a tensor component must be added. For comparison, 
the A-CDM concordance model is also plotted. The middle 
panel shows the expected polarization signal of the best-fitting 
model. Including the tensor components produces B modes in 
the polarization spectra. In the bottom panel it can be seen 
that the best-fitting model to the BOOMERanG-98 data is in 
contradiction with the matter power spectrum measured by the 
PSCz survey. 



One can further combine the CMB constraints with 
those obtained from the lum inosity-distance measure 
mcnts of high-z sup ernovae ( |Pcrlmutter et al. 1998 
Schmidt ct al. 199S ) as we do in the bottom panel of 



Figure 1. The inclusion of the SNIa data does not greatly 
affect the constraints on £7a but, as expected, it rules out 
any models with low values of fib- We find that fif, > 0.13 
at the 95% C.L. and an increasing inconsistency of BDM 
models with a flat universe. 

Also plotted in the bottom panel of Figure 1 are the con- 
straints obtained including a p rior Vl^h 2 = 0.019 ± 0.002 
(Buries, Nollett & Turner 2000) inferred by measurements 



of primordial elements and assuming standard BBN. The 
results of the 2 priors are incompatible at more than 3er. 

It is interesting to further explore the compatibility of 
our results with independent observations by considering 
likelihood contours in the (fif,, h) plane. Figure 3 plots the 
results of our analysis with no priors, the results including 
the SNIa data and the results 
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FlG. 3. — The likelihood contours in the (f2s, /i), with the 
remaining parameters taking their best-fitting values to the 
BOOMERanG-98 data. The contours correspond to 0.32, 0.05 
and 0.01 of the peak value of the likelihood, which are the 68%, 
95% and 99% C.L. respectively. The filled contours correspond 
to an analysis with a t > 10 Gyr prior. The dashed, dotted and 
solid contours to the right of the plot are the result when the 
SNIa data is included. The bold black line to the left of the plot 
is the 95% C.L. region from BBN. with the BBN constraint. 

We can see from this figure that purely baryonic models 
are in quite good agreement with BOOMERanG-98 and 
with the BBN constraints. However, when the SNIa data 
is included we are restricted to models with values of Qb 
that are too high and values of h that are too low to be 
consistent at < 3er with both the BBN constraint and th e 
recent HST result h = 0.72 ± 0.08 ( [Freedman et al. 2000] ). 

Quite independently of the prior assumed, we find that 
the age of the universe in this scenario is constrained to be 
to = 24±2 Gyr, nearly double that expected from standard 
CDM. The age of the universe must exceed the ages of the 
oldest globular clusters tec = 14 ± 2 Gyr by an amount 
At ~ 0.5 — 2 Gyr, so the age of the universe required by 
these models seems to be too high unless globular cluster 
formation is delayed for a very long time. 

Collisional damping erases fluctuations in the matter 
power spectrum on scales < l00h~ 1 Mpc. The predicted 
matter power spectrum from the best-fitting model (see 
Figure 2, bottom panel) is in contradiction with the decor- 
related linear po wer spectrum extracted from the recent 
PSCz catalogue ( [Hamilton fc Tegmark 2000[ ). Allowing 
for a shift in the overall amplitude with a bias factor b 
where P{k)pscz — b 2 P(k), we obtain a best-fitting bias 
b = 4.5 with a x 2 — 155 which is an extremely poor fit 
to 22 data points with one free parameter. The theoreti- 
cal variance in baryonic matter fluctuations over a sphere 
of size 8/i _1 Mpc, for example, is erg ~ 10~ 2 to be com- 
pared with the observ ational value <7g = 0.56f2~ 047 ~ 1.4 
( Viana fc Liddlc 1999| ). So, rescaling the galaxy data with 
a massive bias parameter may not be physically meaning- 
ful since no galaxies would have formed in order to allow 
there to be a predicted galaxy power spectrum. 

Various phenomenological mechanisms can be proposed 
in an attempt to resolve the discrepancy with the mat- 
ter power spectrum and render BDM models viable. One 
such possibility is to increase the spectral index of the 
primordial density fluctuations, n s . However, the 'blue' 
models that are able to produce a$ ~ 0.5 require n s > 2 



tion ( {Steinhardt 1995 ). Also, such a high tilt will cause the 
CMB spectrum to become incompatible with observations. 

Another option that could increase small scale den- 
sity perturbations whilst leaving the CMB spectrum un- 
changed is to have a spectral index that varies with scale, 
dn s /dlnfc =fi 0, as predicte d by some inflationary models 
( Kosowsky fc Turner 1995 ). Unfortunately, the variation 



needed, 



1 , is much greater than that predicted by the 



most viable models. 

More recently, there has been some investigation into 
features in t he primordial power spectrum g enerated dur- 
ing inflation (Griffiths, Silk & Zaroubi 2000; ; Knebe, Islam 



& Silk 2001). It is possible that such features could in- 



crease the perturbations on small scales without affecting 
the CMB spectrum, although it is doubtful that a single 
feature would be able to resolve the observed discrepancy 
and such a mechanism is unlikely to be preferable to the 
assumption of CDM itself. 

4. DISCUSSION 

We have examined CMB anisotropies and large-scale 
structure observations in a purely baryonic dark matter 
universe. Our results suggest that a BDM adiabatic uni- 
verse with a power-law primordial power spectrum can 
only reproduce the sub-degree CMB measurements if the 
universe is closed. Although the premise that fi = 1 is 
considered to be one of the basic predictions of the infla- 
tionary scenario, it is possible to construct infl ationary see 



Linde 1995 



which is outside the range 0.7 



1.2 allowed by infla- 



narios that predic t closed universes (see e.g 
Starobinsky 199£ ). It is interesting to note, however, that 
removing CDM immediately forces us to construct a more 
elaborate inflationary model, or, in other words, the most 
simple model needs non-baryonic matter to work. 

When flatness was assumed, we found that ui c d m > 0.04 
at 2cr with L(uj cdm = Q)/L{oj cdm = 0.1) — 10~ 3 . There- 
fore, the flatness constraint with BBN renders ABDM 
models less consistent with the CMB observations than 
the standard ACDM model. Also, ABDM models that 
fit the BBN constraints are degenerate with a non-physical 
region of the parameter space (with very high baryon con- 
tent and small h). When any cosmological prior other 
than BBN is assumed in the analysis, the unphysical re- 
gion turns out to be preferred. Including the SNIa con- 
straints (or f2 m = f2f, > 0.25), for example, makes the 
parameter space incompatible at > 95% confidence with 
both the HST constraint h = 0.72 ± 0.08 and with BBN. 

As expected, all the BDM adiabatic models fail to repro- 
duce both the CMB observations and the observed amount 
of galaxy clustering on 8/i _1 Mpc scales. In order to solve 
all the discrepancies, that are not present in the standard 
ACDM scenario, we need to introduce ad hoc mechanisms 
which are unlikely to be preferable to the assumption of 
CDM itself. Nonetheless, BDM models leave a set of char- 
acteristic imprints such as large scale B-mode polarization 
and no third acoustic peak in the anisotropy spectrum that 
will allow future experiments to further scrutinize this hy- 
pothesis. 
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